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Estimate*  of  Vertical  Eddy  Diffusion  Due 
to  Turbulent  Layer*  in  the  Stratosphere 


1.  INTRODUCTION 

The  present  concern  about  the  vulnerability  of  stratospheric  ozone  to  pollution 
raises  the  question  of  how  material  moves  vertically  in  that  stable  environment. 

In  the  case  of  an  homo'geneouu  fluid  (with  neutral  stability),  turbulence  playB  a domi- 
nant role  in  mixing  and  vertical  transport.  The  situation  is  markedly  different  in 
the  case  of  a stably  stratified  fluid  which,  of  course,  is  the  case  of  interest  in 
studies  of  the  stratosphere.  In  this  case,  the  turbulence  confines  itself  to  thin 
horizontal  layers  separated  by  non-turbulent  layers.  The  turbulent  layers,  which 
are  caused  by  the  Kelvin- Helmholtz  <KH)  instability  and  which  have  the  geometry 
(roughly)  of  a flat  disc  ("blini"  or  pancake),  usually  occupy  a Bmall  fraction  of  the 
fluid  volume  (for  example,  5 percent)  and  in  the  case  of  the  stratosphere  their 
vertical  extent  is  presumably  of  order  200  meters.  In  contrast,  their  horizontal 
extent  is  typically  in  the  range  of  a few  tens  of  kilometers.  The  same  sort  of  situa- 
tion prevails  in  the  upper  ocean  but  there  the  vertical  dimension  is  reduced  to 
order  10  centimeters. 

The  vertical  eddy  diffusion  coefficient,  K , for  the  stratosphere  over  long  time 
averages  (based  on  radioactive  fal'out  measurements)  ranges  from  about  0, 1 m*/sec 
at  the  equator  to  1.  0 m2/sec  at  t'  ; poles.  This  estimate  (Junge), 1 of  course, 

(Received  for  publication  1 February  197  9) 

1.  Junge,  C.  E,  (1963)  Air  Chemistry  and  Radioactivity.  Academic  Press, 

New  York,  p 260, 
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includes  not  only  the  effects  of  the  above-mentioned  small  scale  turbulence,  but  it 
also  includes  global  circulation  and  other  effects  as  well.  An  important  but  un- 
answered question  is  whether  or  not  the  small  scale  turbulence  plays  a significant 
role  in  the  overall  value  of  K . 

One  method  to  estimate  K due  to  turbulence  depends  upon  estimates  of  the 
dissipation  rate,  e , basod  in  turn  upon  the  spectral  analysis  of  in  situ  velocity 
fluctuations.  Lilly2  has  done  this  and  found  K0~O.  01  m‘/sec  which  is  very  small. 
The  purpose  of  the  present  report  is  to  closely  examine  a newer  method  to  estimate 
K (due  to  small  scale  turbulence)  which  is  based  on  measurements  of  the  vertical 

v 

profile  of  the  horizontal  mean  winds  made  in  conjunction  with  measurements  of  the 
vertical  temperature  profile.  Suoh  wind  measurements  have  been  made  in  our  field 
experiments  by  means  of  rocket  trails  and  by  others  (Van  Zandt  et  al)3 4 5  using  radar 
techniques.  Temperature  profiles  are  obtained  from  measurements  from  free 
balloons  carried  out  simultaneously.  Rosenberg  and  Dewan*  were  the  firBt  to  des- 
cribe this  technique  and  the  reader  will  find  in  that  report  an  extensive  review  of 
the  literature  leading  to  the  main  concepts  involved. 

More  precisely,  the  goal  of  this  report  is  to  make  explicit  all  the  assumptions 

behind  our  model  and  to  examine  them  in  detail.  Also,  generalizations  will  be  dea- 

o 

orlbed  which  will  show  when  it  is  permitted  to  use  the  relation  employed  by  Lilly 
(and  in  the  context  of  the  ocean  by  Woods  and  Wiley  ) given  by 


K ■*  P*Kt 
e u 


(1) 


where  P*  is  the  fraction  of  the  vertical  dimension  which  is  occupied  by  turbulent 
flow  and  is  the  effective  average  eddy  diffusion  within  the  turbulent  layers.  It 
will  be  shown  that  when  is  large,  Eq.  (1)  is  not  valid  and  one  must  use 
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IB 
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where  L is  related  to  turbulent  layer  thickness  and  Atj  is  the  average  tlmo  between 
changes  of  the  vertical  Richardson  number  profile.  Eqs.  (1)  and  (2)  will  be  derived 


2.  Lilly,  W. , Waco,  L. , and  Adelfang,  S,  (1975)  Stratospheric  mixing  estimated 

from  hlgh-altltude  turbulence  measurements  by  using  energy  budget  techniques. 
The  Natural  Stratosphere  of  1974,  CIAP  MONOGRAPH  1,  DOT-TST-7  5-51 

pp'6-81  to'B-BO. 

3.  Van  Zandt,  T,  E. , Green,  J.  L. , Gage,  K,  S. , and  Clark,  W.  L.  (1078)  Vertical 

profiles  of  refractlvity  turbulence  structure  constant,  Radio  Scl,  JiShB  19-829. 

4.  Rosenberg,  N.  W. , and  Dewan,  E.  M,  (1975)  Stratosphere  Turbulence  and 

Vertical  Effective  Diffusion  Coefficients,  Al‘'tJRL-TR-75-ffB19,  ErP  No.  535. 

5.  Woods,  J.D.,  and  Wiley,  R.  L.  (1972)  Billow  turbulence  and  ocean  microstructure. 

Deep  Sea  Research  and  Oceanic  Abst.  19s 87 -121. 


from  two  approximations  resulting  from  the  vertical  stack  model  which  was  first 

4 

described  in  ltosenberg  and  Dewan. 

In  addition  to  discussing  certain  subtleties  of  the  vertical  stack  model  us  well 
some  details  of  data  analysis,  this  report  compares  the  behavior  of  turbulent 
layers  in  the  atmosphere  to  those  in  the  upper  ocean.  The  latter  Is  of  interest 
because  the  upper  ocean  is  much  more  accessible  than  the  stratosphere  and  acts, 
to  some  extent,  as  a "laboratory  model"  of  it.  Finally,  a list  of  Important  experi- 
mental questions  raised  by  the  model  is  given.  Preliminary  results  using  the 

A 

present  approach  were  given  in  Rosenberg  and  Dewan  hereafter  designated  by  I; 

o 

and  an  estimate  of  K ~0,  3 m / sec  was  given.  If  such  an  estimate  were  to  be 
validated  by  future  measurements,  then  turbulence  would  be  important  in  vertical 
transport  in  the  stratosphere.  At  present,  however,  the  question  remains  open. 


2.  REMARKS  ON  THE  SIMILARITIES  ANI)  DIFFERENCES  IN  THE 
TURBULENCE  OF  THE  UPPER  OCEAN  AND  STRATOSPHERE 

a K 

The  studies  by  Woods  and  Woods  and  Wiley  of  the  thermocllne  and  layered 
turbulence  to  be  found  there  are,  as  has  been  mentioned,  very  relevant  to  stratos- 
pheric research  (cf.  I).  The  thermocllne,  or  region  where  there  is  a stable  tempera- 
ture gradient  near  the  ocean  surface,  seems  to  be  divided  Into  regions  which  are 
called  "sheets"  across  which  are  very  sharp  density  gradients,  separated  by  "layers" 
where  the  gradients  are  much  smaller.  Turbulence  occurs  In  the  vicinity  of  the 
sheets  whereas  the  layers  can  be  considered  to  be  effectively  laminar. 

An  important  role  seems  to  he  played  in  thermocllne  turbulence  by  "buoyancy 
waves,"  sometimes  called  internal  waves.  TheBe  waves  (as  was  shown  theoretically 
by  Phillips)  can  enhance  local  mean  shears  In  a manner  which  results  In  lower 
R,.  When  R,  < 0.  26,  turbulent  breakdown  can  take  place,  and  when  turbulence  coni- 
nietices,  as  described  by  Woods  and  Wiley,  an  hysteresis  phenomenon  Is  set  into 
motion,  The  turbulence  causes  the  layer  to  thicken,  but  this  thickening  causes  Rj 
to  increase.  Turbulent  entrainment  is  the  cause  of  the  thickening  and  the  latter 
can,  on  the  basts  of  energy  considerations,  continue  until  R^  = 1.  Once  Rj  » 1, 
however,  the  turbulence  can  no  longer  extract,  energy  from  the  mean  shear  and 
hence  R(  = 1 represents  the  "cut  off.  " Thus,  during  the  "active  life"  of  the  turbu- 
lence, the  layer  expands,  and  at  ■*  1 motion  decays  and  this  decay  endures  over 
a period  of  time.  In  the  case  of  the  ocean,  the  layer  expands  by  a factor  of  four 
(during  the  active  life).  This  was  determined  by  observation  by  Woods  and  Wiley8 

0,  Woods,  ,1.  D.  (1968)  Wave-induced  shear  Instability  lr.  the  summer  therinocllne, 

J.  Fluid  Mech,  ^7  91-800, 

7,  Phillips,  O.  M.  (1989)  The  Dynamics  of  the  Upper  Ocean  (Cambridge  tJnlv.  Press). 
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who  also  estimated  It  on  the  basis  of  a simple  model  which  assumes  that  the  hori- 
zontal velocity  difference  across  the  layer,  W,  as  well  as  the  temperature  differ- 
ence, ^O,  remain  constant  during  expansion.  After  layer  expansion,  a "collapse" 
takes  place.  This  presumably  is  due  to  the  effect  of  incomplete  mixing  followed  by 
the  return  of  fluid  parcels  to  their  stable  depthB  (Koop). 8 

The  appearance  of  tropospheric  clear  air  turbulence  (CAT)  on  radar  has  a 
striking  similarity  to  that  seen  in  the  ocean  and  laboratory.  But  there  are  also 
important  differences.  For  example,  in  the  case  of  Browning  and  Watkins'  obser- 

D 

vation,  there  is  a complete  absence  of  collapse  subsequent  to  the  expansion.  The 
expansion  effects  themselves  don't  seem  to  be  visible  on  the  radar.  In  addition,  as 
will  be  shown  below,  one  would  expect  that  in  the  case  of  the  atmosphere  there  is 
much  less  expansion  than  in  the  ocean,  and  from  the  above  remarks,  it  should  be 
clear  that  this  in  turn  would  help  explain  the  lack  of  collapse.  Specifically,  with  less 
turbulent  entrainment,  there  would  be  less  fluid  to  mix  and  thus  mixing  could  be 
more  thorough,  and  hence  there  would  be  far  more  "homogenization"  and  less  migra- 
tion of  parcels  to  their  stable  altitude. 

In  order  to  explain  the  hypothetical  lack  of  turbulent  layer  spread  in  the  atmos- 
phere, we  note  that  the  above-mentioned  assumption  of  constant  & 6 over  time  across 
the  layer  made  in  the  case  of  the  ocean  is  probably  invalid  for  the  case  of  atmos- 
pheric  turbulence.  In  the  case  of  the  thermocline  the  temperature  profile  is  deter- 
mined by  physical  mixing  of  water  (Phillips)7  and  therefore,  between  mixing  events 
it  remains  unchanged.  In  contrast,  the  6 profile  for  air  (here  we  Identify  8 with 
potential  temperature)  is  due  predominantly  to  radiative  effects  and  between  the 
mixing  events  the  "steps"  in  the  profile  would  tend  to  become  smoothed  out.  It  Is, 
therefore,  more  appropriate  to  roplace  the  constant  A0  condition  with  a constant 
d0  /dz  or  vertical  gradient  (or  buoyancy  frequency,  N0)  condition.  When  this  is 
taken  into  account  and  when  we  furthermore  replace  the  R.  » l cut  off  value  with  the 
one  advocated  by  Garrett  and  Munk  (which  was  due  to  Thorpe)  namely  R^  ■ 0,4, 
then  It  can  be  shown  (see  I)  that  the  expansion  factor  1b  a mere  1.  26  in  contrast  to  4 
(R^  ■ 1 would  have  given  an  expansion  factor  of  2). 

Another  point  raised  by  Woods  and  Wiley5  is  that  "sheet-ensembles"  might  be 
generated  by  repeated  billow  events.  This  in  turn  raises  the  possibility  that  the 
same  sort  of  thing  occurs  in  the  stratosphere. 


0.  Koop,  C.G.  (19? 8)  Instability  and  Turbulence  in  a Stratified  Shear  Layer, 

USCAE  134,  Univ.  of  Southern  California,  School  of  Engineering,  Dept.  of 
Aerospace  Engineering. 

9.  Browning,  K.A.  (1971)  Structure  of  the  atmosphere  in  the  vicinity  of  large 

amplitude  Kelvin -Helmholtz  billows,  Roy.  Met.  Soc.  Quart.  J.  J37 : 263-290. 

10.  Garrett,  C.  , and  Munk,  W,  (19  ? 3)  Oceanic  mixing  by  breaking  internal  waves, 
Deep  Sea  Research  JJgj B 2 3 - B32, 
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3.  THE  VERTICAL  STACK  MODEL  FOR  Ke 

3.1  Assumptions 

First,  we  shall  assume  that  there  is  no  vertical  transport  taking  place  between 
the  turbulent  layers.  Certainly  a significant  amount  of  transport  by  molecular 
diffusion  (as  compared  to  turbulent  diffusion)  can  be  ruled  out  on  the  basis  that  the 

molecular  diffusion  is  of  order  10“^  m^/sec  in  the  lower  stratosphere  whereas  K 

2-22  -12  e 
(according  to  Lilly)  is  10  m /sec,  or  (according  to  I)  10  m /second.  The 

second  assumption  we  make  is  that  turbulence  occurs  randomly  with  altitude  and 
time,  With  no  information  to  the  contrary,  this  is  a valid  a priori  assumption. 

Next  we  assume  that,  before  the  turbulence  decays  in  the  layer,  total  mixing  takes 
place.  This  assumption  receives  support  from  two  entirely  independent  pieces  of 
experimental  information,  The  first  is  described  in  the  report  by  Mantis  and  Pepin11 
where  it  was  shown  that  in  the  stratosphere  there  are  layers  where  the  temperature 
gradient  with  respect  to  altitude  is  nearly  equal  to  the  adiabatic  lapse  rate.  This 
would,  of  course,  be  consistent  with  the  idea  that  total  mixing  occurs  because  such 
mixing  would  presumably  leave  behind  an  adiabatic  lapse  rate,  The  second  piece 
of  evidence  is  from  Browning  and  Watkins  and  also  from  Atlas  et  al.  These 
authors  looked  at  turbulent  layers  ("CAT")  in  the  troposphere  with  radar  and  they 
found  that  after  the  turbulent  activity  had  presumably  decayed  (following  the  In- 
stability) there  remained  two  widely  separated  layers  of  very  strong  reflection  each 
located  at  an  edge  of  the  original  layer.  This  Indicates  a strong  gradient  of  tempera- 
ture at  each  edge  which  is  precisely  what  would  be  found  if  essentially  total  mixing 
occurred  within  the  layer.  In  a sequel  to  the  present  report,  1 shall  give  a theoretical 
argument  which  gives  still  further  justification  to  the  total  mixing  assumption.  (The 
case  of  small  mixing  will  also  be  examined  below  as  already  mentioned.  ) 

Finally,  it  is  assumed  that  the  horizontal  rearrangement  of  the  layer  and  trans- 
ported material  will  huve  little  effect  on  vertical  transport.  This  last  assumption 
is  not  unusuul  in  oceanography,  und  it  permits  the  use  of  a simple  one-dimensional 
model  (see  Garrett  und  Munk).  10 

3.2  The  Model  for  Kv 

The  elusalcul  eddy  diffusion  concept  was  formulated  for  the  case  where  tur- 
bulence Is  relatively  uniform,  That  is  to  say,  it  would  not  apply  to  our  case  where 


11.  Mantis,  M.T.,  und  Pepin,  T.J.  < 107  1)  Vertical  temperature  structure  of  the 

free  atmosphere  at  mesoscale,  J,  Cleophyo,  lies.  20; 8621- BO 2(1. 

12.  Atlas,  II.,  Metcalf,  J.,  Richter,  .].,  and  Gossard,  IS.  (197  0)  Tile  birth  of 

"CAT"  and  microscale  turbulence,  J.  of  the  Atm.  Set.  903- 913. 
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certain  regions  of  the  fluid  are  turbulent  while  others  are  laminar.  The  usual 
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approach  (see  for  example,  Pasqulll)  is  briefly  as  follows. 

Let  C be  the  concentration  of  material  as  a function  of  position  and  consider 
the  vertical  flux,  Q,  due  to  the  turbulence.  Letting  w'  be  the  fluctuation  of  vertical 
velocity  andc'  the  fluctuation  in  concentration,  then  Q = w'  Cf  where  the  overbar 
signifies  a suitable  average.  C'  is  given  by  C'  * -£8  C/8  2 where  £ was  originally 
called  a "mixing  length"  but  is  now  considered  to  he  a characteristic  spatial  dimen- 
sion of  the  turbulence.  Kg  is  then  defined  by  analogy  with  Fourier's  heat  transfer 
equation 

Qi  -K  8 TT/8  Z (3) 

and  thus  Ke  ■ w'£. 

This  approach  fails  for  the  case  we  are  considering,  therefore  we  re-examine 
it  in  detail  in  order  to  arrive  at  an  extension  of  the  eddy  diffusion  concept  which 
can  be  of  use  in  stratified  fluids, 

In  Figure  1 the  fluid  is  modeled  by  a medium  located  between  a source  at  the 
top  and  a sink  at  the  bottom,  Within  this  medium  are  few  randomly  spaced  mixing 
layers  and  we  assume  that  the  Initial  concentration  profile  is  linear  as  indicated 
in  Figure  la.  Figure  lb  indicates  the  effect  upon  the  profile  of  the  mixing  layers 
where  the  assumption  of  total  mixing  aocounts  for  the  vanishing  gradients  within 
the  layers.  Note  that  there  results  a net  downward  transport  from  this  mixing, 

This  fact  is  essential  for  the  understanding  of  this  process  of  vertical  transport. 

Note  also  that  if  the  layers  were  to  be  locked  to  this  initial  configuration,  no  further 
vertical  transport  through  the  medium  would  be  possible,  We  shall  assume  that 
whenever  a layer  is  in  contact  with  either  the  source  or  the  sink,  the  mixing  action 
will  convert  its  concentration  to  that  of  the  source  or  to  zero  in  the  case  of  the  sink. 
In  the  context  of  stratospheric  transport  the  "instantaneous  sink"  is  justified  from 
the  fact  that  the  residence  time  is  of  the  order  of  years  for  the  stratosphere,  where- 
as in  the  troposphere  It  is  given  in  terms  of  weeks,  Storm  systems,  convection, 
and  other  transport  mechanisms  operate  In  the  latter  and  these  are  responsible  for 
the  great  disparity  in  residence  times. 

Next  we  imugine  that  the  medium  in  Figure  l is  subjected  to  a succession  of 
profiles  of  randomly  spaced  mixing  layers  of  variable  thickness,  Figure  3 shows 
a computer  simulation  of  this  process  after  a significant  number  of  profiles.  As 
cun  be  seen,  the  C (2)  profile  departs  from  the  original  UneHr  shape  and  becomes 
very  choppy.  Thus  this  form  of  transport  is  very  "lumpy"  and  the  "diffusion" 
represented  by  the  Kq  derived  below  must  be  viewed  in  this  context. 


13,  Pasqulll,  F.  (1002)  Atmospheric  Uiffusion,  Van  Nostrand  Co,  Ltd, 
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Figure  la,  Vertical  Stack 
Model  of  Turbulent  Mixing 
Layers.  Ordinate  is  altitude, 
slanted  line  is  concentration 
profile  and  abscissa  desig- 
nates amount  of  concentration 


Figure  lb.  Effect  of  Mixing 
Layers  on  Concentration 
Profile.  Total  mixing  inside 
layers  is  assumed.  Net 
downward  transport  from 
source  to  sink  is  shown 


CONCKNTNATIONlCu) 

Figure  2.  Computer  Simulation  of  the  Effect  Upon 
Linear  Concentration  Profile  (dotted  line)  by  a Series 
of  "Time  Frames"  of  Turbulent  Mixing  Layer  Profiles. 

In  this  simulation,  statistics  of  thickness  and  occurrence 
are  based  on  experimental  wind  profiles  in  the  stratos- 
phere. Discrete  nature  and  step  like  structure  is 
illustrated 
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Presumably,  in  the  actual  stratosphere,  the  mixing  layers  occur  at  random 
altitudes,  with  random  thickness  at  random  times.  To  model  this  in  a simple 
manner  we  conceptualize  the  process  as  occurring  in  discrete,  evenly  spaced  time 
frames  with  a single  profile  of  mixing  layers  in  a given  frame.  The  duration  of 
these  time  frames  we  denote  by  At^. 

From  Eq.  (3)  we  define 

Ke  a -Q/(0T?/8Z)  . (4) 

We  must  now  determine  the  value  of  the  flux  Q.  This  flux  is  extremely  sporadic 
in  the  sense  that  no  material  is  "dumped"  out  of  the  bottom  of  the  stack  (Figure  1) 
until  a layer  forms  there.  We  thus  envision  the  process  in  a manner  which  regards 
it  only  over  a long  period  of  time  and  in  an  averaged  out  sense.  On  average,  the 
stack  is  in  "steady-state"  with  as  much  material  entering  the  top  as  leaving  the 
bottom.  The  profile  C (Z)  is,  on  average,  a roughly  linear  profile. 

To  calculate  Q we  define  Atb  as  the  average  time  between  the  occurrences  of 
layers  at  the  bottom  of  the  stack  (in  contact  with  the  sink),  This  can  be  related  to 
Atj  from  the  definition  of  P>'\  Thus 

P*  ■=  Atf/Atb  . (5) 

Since  P#  is  very  small  (for  example,  0.  05),  Atb  In  general  will  be  much  larger  than 
Atf.  Next  we  define  the  symbolAn  to  equal  the  amount  (on  average)  of  the  material 
in  the  bottom  layer  whioh  Is  then  "dumped  out"  into  the  sink  due  to  mixing.  Let  L 
in  Figure  3 represent  the  average  thickness  of  the  bottom  layer.  Assuming  the 
linearity  of  C (Z)  the  average  concentration  in  the  layer  Is  1/2  • L • 0 C/B  Z and  from 
this  one  obtains 

An  - (A  L)  ’ (L/2)  0 C/8  Z (6) 

where  A is  the  area  of  bottom  surface  and  A L Is  the  volume  of  matter.  Q there- 
fore is  now  obtained  from  Eq.  (6)  as 


(the  minus  occurring  because  the  flow  is  downwards). 
Therefore,  from  Eq.  (4) 


,\ouRce 


/'  TRANSFER  OF 

MATERIAL  FROM  TURBULENT 
LAYER (L)  TO  SINK 

An*  NO  OF  Kg  TO  SINK*  L -Cay -A 
Coy  ■ dc/dj  ■ L /i 

C0*Cma»  OF  BOTTOM  LAVER 


Figure  3.  Diagram  for  Ke 
Vertical  Stack  Model  Relating 
Wind  Profile  Statistics,  f>  (Llto 
Vertical  Transport.  Randomness 
and  stationarity  relate  the  dumping 
effects  of  the  bottom  layer  to  the 
entire  stack  structure 


i mvr  L I <rrr> 


As  mentioned,  the  value  of  Lr  must  relate  to  data  by  means  of  some  sort  of 
average;  however,  before  we  see  how  this  is  carried  out  in  practice  it  La  Important 
to  understand  the  subtle  meaning  of  L.  For  example,  it  is  definitely  not  to  be  under- 
stood as  the  average  thickness  of  a turbulent  layer  in  the  fluid.  Instead  it  is  the 
average  thickness  of  a turbulent  layer  whioh  extends  above  a previously  given 
altitude.  ''  Recall  that  the  bottom  of  the  stack  is  at  a given  fixed  location,  The 
next  section  will  reflect  this  consideration  by  the  way  in  which  the  data  is  processed. 


3.3  L2  and  Analyili  of  Data 

— 2 

We  now  consider  the  experimental  determination  of  h to  be  used  in  Eq.  (8)  to 
obtain  Ke«  The  primary  data  consists  of  vertical  profiles  of  the  mean  horizontal 
winds.  In  addition  one  also  needs  the  temperature  profiles.  The  data  are  combined 
to  give  the  HichardBon  number  (Rj)  profiles  where 

u . (g  a $ im  /fdui2 


and  g is  the  acceleration  of  gravity,  9 the  potential  temperature,  8 its  average, 
end  “ is  the  average  horizontal  wind  velocity.  In  order  to  determine  the  potential 
existence  of  turbulence,  the  usual  criterion  of  R^  < 0.  2b  is  used  (cf.  Ref.  I). 

The  purpose  of  the  data  analysis  is  to  search  the  wind  profiles  for  supercritical 
shears  (R^  < 0.  25).  To  illustrate  the  technique  we  describe  in  detail  what  was  done 


Thus  on  average  the  layer  thickness  is  roughly  twice  the  average  value  of  L. 
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to  obtain  the  results  given  in  Ref.  I.  In  that  case,  data  were  available  at  29  m 
altitude  Intervals.  Since  no  temperature  data  were  available,  we  used  a standard 
atmosphere  and  set  the  numerator  of  Eq,  (9)  to  that  (constant)  value.  A computer 
program  searched  the  wind  profiles  for  unstable  zones  as  follows.  For  a given 
profile  the  search  began  at  the  bottom  point.  The  shear  between  this  point  and  the 
one  a step  higher  was  calculated.  If  Rj  > 0,  25,  the  sample  was  counted  as  non- 
orltloal  (non-turbulent)  and  the  program  considered  the  second  point  up  ae  its  initial 
altitude  and  repeated  the  measurement  as  before.  When  an  altitude  point  was 
reached  where  the  shear  between  that  altitude  and  the  subsequent  one  was  super- 
critical, the  program  would  then  calculate  the  shear  between  the  Initial  altitude  and 
one  two  steps  higher  (50  m),  then  three  steps  (75  m)  and  on  until  an  altitude  span 
was  reached  where  R^  > 0.  25.  The  program  would  then  return  downwards  one  step 
and  record  the  altitude  difference  as  a value  of  L.  When  all  these  values  were 
obtained  they  were  counted  up  in  the  form  of  a cumulative  fraction,  Pj  (L), 

Pj(L)  is  defined  as  the  fraction  of  altitude  samples  which  had  a supercritical 
shear  between  the  given  altitude  and  one  a total  distance  L above  it  (or  greater). 

In  other  words 

L max 

PT(L)"E  AP(L)  (10) 

1 L 

where  A P (L)  represents  the  fraction  of  samples  with  L between  the  value  of  L and 
L + AL,  where  A L " the  altitude  interval  between  samples  (for  example,  20  m). 
From  Eq.  (10) 

L max  L-AL 

Pv(L)»  E AP  (L)  - E AP  (L)  (11) 

1,  min  L<  min 

where  L min  = AL  and  L max  is  the  maximum  span  of  super- critical  shear.  Re- 
calling that  P*>  is  the  average  fraction  of  vertical  profiles  that  are  occupied  by 
unstable  shears  (potentially  turbulent)  we  have 

L max 

P*  * E AP  (L>  . (12) 

L min 

From  Eqs.  (11)  and  (12) 

L-AL 

P.  <L)  t=P*-  E AP  (L).  (13) 

1 L min 

(L  can  be  replaced  by  mA  L whore  m is  an  integer  from  1 to  some  upper  limit. ) 

AP  (L)  is  then  calculated  from  Pj  (L)  by 
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AP(mA  L)  = Pr  (rrsAL)  - P [(m  -c  DAL]  (14) 

as  can  be  seen  from  Eq.  (13). 

— 5 

Next,  one  needs  a probability  in  order  to  calculate  L and  for  this  we  define 

P (L)  ■— jffil  (15) 

which  represents  the  fraction  of  potentially  turbulent  zones  found  which  have  the 
span  L.  Finally 

— w m ■ max  n „ 

L*  ■ T.  (mAL)  P(mAL).  (1.6) 

m » 1 


K is  then  determined  from  Eq.  (8)  or,  using  Eq.  (B) 
0 


K 


e 


P*  L‘ 

Tzrtj- 


(17) 


Ref.  I carried  out  these  procedures,  Figuro  4 shows  P (L)  based  on  10, 000 
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samples  of  wind  profile  data  published  by  Miller,  Henry  et  al  at  NASA.  To  sum- 

*TT  4 2 

marizc  the  findings  of  Ref.  I,  we  used  Eq.  (17)  to  obtain  L » 1,  2 X 10  m , 

Eq.  (12)  to  obtain  P’:'  * 4,  8 X 10’2,  and  Eq.  (17)  then  gave  Ke  - 0,  10  m2/seoi  (Atf«= 
1500  sec  was  used  and  this  will  be  discussed  at  length  below).  The  previously 
discussed  expansion  factor  (due  to  layer  spread)  results  in  a factor  1,6  increase  in 
K and  thus  we  arrived  at  the  estimate  of  K„  ~ 0,  3 m2/second. 

0 R 


3.4  Rxtcnilon  to  the  ()«■«  of  Smell  Mixing 

Having  obtained  Eq.  (18)  on  the  assumption  that  complete  mixing  takes  place 
within  a turbulent  layer  we  now  turn  to  the  opposite  case  where  there  is  only  a small 
amount  of  mixing  and  derive  Eq.  (1)  which  was  used  by  Lilly  et  ul,  from  the  verti- 
cal stack  model,  is  defined  as  the  eddy  diffusion  coefficient  within  a mixing 
layer.  Turning  again  to  Figure  3 and  the  bottom  layer  of  the  stack,  we  estimate  An 
in  the  new  context,  Let  Atm  represent  the  duration  of  the  mixing  and  let  Kj  bo  re- 
garded as  constant  during  this  time  and  zero  ufterwards,  (This  simplification  will 
be  made  less  artificial  below  by  considering  the  case  of  monotonically  decreasing 
in  time  without  limit  for  actual  mixing  time, ) Let  nQ  be  the  amount  of  material 
inside  the  layer  locuted  over  a unit  horizontal  area  at  t ■ 0 and  let  n^t  be  this 
quantity  at  t = At^.  Then,  ^ 


14.  Miller,  R,  W,,  Henry,  K.  M,  , and  Rowe,  M.CL  (1965)  Wind  Velocity  Profiles 

Measured  by  the  Smoke-Trail  Method  ut  Wallops  Island.  Vlrelnln  ( l&BO- 1982) 
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Thus,  to  estimate  An  one  must  know  n as  a function  of  time  as  it  diffuses  out  of 
the  layer  n.  The  definition  of  n is  the  amoui.'  of  material  inside  the  layer  at  time  t 
over  a unit  horizontal  area.  As  seen  in  Figure  3,  C0  is  the  maximum  concentration 
of  the  layer  thus 


C (2) 


(10) 


The  heat  equation  which  can  be  derived  from  Eq.  (4)  with  KQ  replaced  by  KL  is 


8 C 


t r 


(20) 


To  complete  the  boundary  value  problem,  it  la  assumed  that  the  top  of  the  turbulent 
layer  Is  "non-conducting"  that  is,  the  normal  derivative  of  C with  respect  to  2 at 
the  top  surface  will  be  zero  (that  is,  at  Z ■ L)  for  all  times  after  t * 0,  and  that  at 
Z ■ 0,  C in  held  to  the  value  zero  (being  in  direct  contact  with  the  sink). 


0 80  100  800  800  400  m 
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Figure  4,  Histogram  for  Normalized 
Distribution,  F(L) 
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This  sort  of  problem  is  well  known  in  heat  conduction,  and  the  solution  is 
obtained  in  the  classic  manner  of  separation  of  variables  and  Fourier  aerleB  expan- 
sion (see  for  example,  Ingersoll  et  al,  p.  126).  15  One  thus  arrives  at 


C<2,  AtM> 


C S 
0 m«  1 


exp 


m' 


KLAtM 


(21) 


where  b is  the  Fourier  coefficient  given  by 
m 


bm“  7 


/ 2Z  »in  aZ  ^ ^ sinaZdZ 

o ■ 1/9  L 1 J 


where  l ■ 2L  and  a ■ (m  II /0).  Thus 

- sin  <m  It). 


4_ 

Ml’ 


b it  2 *in 

m , m“ 


The  concentration  profile  as  a function  of  time  is  then 

c (*-»  ■ c0  «P  j - (t1)  KL4,m  I 

j"a  ,Ln  ^—3—)  - "in  <«.T1)J  Bin  ^ 

The  value  of  n (t)  can  be  obtained  from 

1/2 

n(t)  ■ / C (Z.t)  dZ  , 
o 


(22) 


(23) 


(24) 


(28) 


Thus 

n(t)  = I (-4)  (1-cos  ^)  (2  sin  - sin  mfl)> 

It3  m-1  ma 

0*P  ^ ” (“7“ ) KLAt!Vlj 
and  at  t * 0 


(26) 


15.  Ingersoll,  L.  R. , Zobel,  O.  J. , and  Ingersoll,  A.  C.  (1954)  Heat  Conductkm. 
University  of  Wisconsin  Press  (Madison,  Wise. ). 
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(27) 


ir ! 
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frifr 


KV.ll 


lr 


■ ‘ i ! 


r1 

[Ur 


vr 


"o  l 


(cf.  Kq.  (7)  where  CQ  ■ (j-|  • L)  and  2L). 

Inserting  Eqs,  (27)  and  (26)  into  Eq,  (18)  and  using  t ■Atj^ 


~K~ 


_ c j r 
n . o 

_ . -p 


1 


i - 1*  r — «xp 
n3  mol  m° 


(28) 


where 

i 


(m)  «J^2  sin  ^ - sin  (m  !1)j  £ 1-cos  . 


(29) 


This  function,  g (m),  takes  on  the  values  2,  0,  -2,  0,  and  so  on,  as  m takes  on  the 
values  1,  2,  3,  4,  etc. 

To  obtain  K in  the  case  of  arbitrary  mixing,  we  insert  Eq.  (28)  into  Eq.  (7) 


and  replace  dC/dZ  by  C /L  tn  Eq,  (4)t 


, 2 f ,R  « 

Ke  ” TSTT  1 ’ ^ ^ . ZT  exp 

b . u m<=  l m 


1 


<<“']■  (3°» 


We  first  specialize  this  most  general  form  to  the  largo  mixing  case  again  by  taking 

(31) 


KLAtM  > * 


(srit) 


Under  this  condition,  Eq.  (30)  becomes  Eq.  (17)  losing  Eq.  (3)]  again,  thuH 
there  is  consistency.  The  case  of  Interest  now  is  the  small  mixing  condition  of 

2 


KLMM<< 


(iA) 


(32) 


It  can  easily  be  shown  that  all  the  non-zero  terms  above  the  first  in  the  sum  given 
in  Eq.  (30)  are  negligible  in  this  problem,  Also  (32/d3)  * 1.0.  Thus  Eq,  (30)  can 
always  be  replaced  by 


(33) 


for  all  practical  purposes,  Now,  setting  m «=  1 in  Eq,  (32)  allows  us  to  expand  the 
exponent 
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2 4 2 

This  we  take  for  the  maximum  permissable  value  of  K.  . Using  the  L -1.21X10  m , 
At..  - 1,  500  sec  (say)  we  arrive  at 


* 2,  1 H tn  /sec 


for  Eq,  (32).  The  values  reported  by  Lilly  et  al  are 


0.  28  < KL  < 1.0  m /sec  . 

Thus  his  results  are,  In  this  sense,  self* consistent.  On  the  other  hand,  the  values 

18  2 

reported  by  Zimmerman  et  al,  namely  40  < < 380  m /sec,  violate  Eq.  (41). 

The  issue  seems  not  to  be  settled,  especially  if  measured  values  of  are  not 
known  to  be  the  "initial  values"  in  all  cases. 

3,5  Growth  Timoi,  Decay  Timer,  and  Vertical  vi  Horizontal 
“Probability  of  Turbulence" 

The  "growth  time"  At.  is  related  to  the  exponential  time  constant  for  the  in- 
* fl 

Stability.  T„>  The  latter  is  given  in  the  form  of  the  following  relation  by  Woods 

- _ 1.  3 (shear)'1 


I | 


which  is  valid  for  0,  05  < ltj  < 0.  20.  The  crucial  parameters  are  Kj  and  shear. 
In  contrast,  the  decay  time  could  be  estimated  from 


AtM-3L 


where  L is  now  assumed  to  be  the  outer  scale  of  the  turbulence,  and  and  are 
the  initial  and  final  fluctuation  velocities.  Eq,  (43)  has  been  derived  from 


d (3/2  u2) 


S' ' I 


where  (3/2  u ) is  the  kinetic  energy  of  the  fluctuations  (u  here  is  to  be  understood 
as  the  root  mean  square  velocity)  and  t is  the  dissipation  rate  taken  equal  to  u / L. 
The  outer  scale  L is  assumed  constant  in  thiB  calculation. 


In  order  to  compare  AtM  with  r , we  must  first  have  values  for  u^  and  uf.  Let 
us  choose  Uj  ■ e"1  Uj.  Thus  the  "enfolding"  mixing  time,  Atj^  is 


18.  Zimmerman,  S,  P, , and  Loving,  N,  (1875)  Turbulent  dissipation  und  diffusivlties 
in  the  stratosphere  using  Htchardson’s  technique.  The  Natural  Stratosphere 
of  1974,  CIAP  MONOGHAPH  1,  Pinal  Report,  DoT'-DT'AP;  DO'TT^’rTr- 5T7 
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We  now  alter  the  form  of  so  that  It  can  be  compared  to  Eq.  (45) , To  do  this 
we  shall  replace  the  shear  in  Eq,  (42)  with  a relation  involving  u1  and  use  an  appro- 
priate value  for  R^.  In  a future  report  concerning  the  degree  of  mixing  subsequent 
to  K-H  breakdown  it  will  be  shown  on  the  basis  of  a simple  energy  argument  that 

/ 1 * Ri  \1/2 

u,  - SL, I — j-i  1 (46) 


where  S is  the  shear  across  the  layer  of  thickness  L. 

Q 

To  choose  R^  we  turn  to  the  observations  of  Browning,  In  view  of  the  fact  that 
L is  often  less  than  his  200  m vertical  distance  resolution  and  in  view  of  his 
range  of  R^  Just  before  the  billow  event  on  Radar  (0,  IB  - 0.  3)  we  arbitrarily  set 
* 0,  15  as  being  a "typical"  value.  Eq,  (4B)  can  be  used  to  obtain  S as  a function 
of  u^  and  L 


,„ui  / Y~ 
b"-n  y}Tl^ 


Putting  this  into  Eq,  (42)  and  again  using  R.  » 0,  IB  we  obtain 

which  should  be  compared  to  Eq.  (45).  Thus  it  is  not  unreasonable  to  expect 


AtM-*tf. 


In  other  words,  despite  a certain  degree  of  arbitrariness  In  the  choice  of  R^.  and  so 
on,  Eq.  (50)  is  certainly  not  unreasonable. 

As  a further  check  on  overall  self-consistency,  let  us  compare  theAtf  obtained 
somewhat  experimentally  in  I with  the  following  estimate  of  (Atf  - 1500  sec). 

Up  to  this  point,  AtM  was  used  together  with  a fixed  value  of  K-L  in  our  calculations. 
We  now  let  be  variable. 

In  this  case,  a good  ''definition"  for  AtM  would  be 


x.v.,.,. aa, 


H,  4*m  ■ I kl 

o 

where  la  the  value  of  at  the  beginning  of  the  turbulence,  Here  y will 
depend  on  the  velocity  fluctuations  and  we  set 


KL(t)  ■ u (t)  L 


t«  3 L 


(m  ’ u“) 


which  Is  derived  in  the  same  manner  as  Kq.  (43).  One  can  solve  for  u (t)  to  obtain 


3L  + u. 


where  ur  is  the  velocity  in  the  first  stage  of  decay, 

* 17  3 

As  pointed  out  in  Tennekes  and  Lumley  however,  £ * u /L  is  no  longer 

appropriate  after  the  Reynold's  number  decays  below  about  10.  At  that  point, 
, _ C ' v U2 


becomes  more  appropriate.  To  match  these  two  expressions  for  £ at  lt0  **  10  we 
must  set  C'  ■ 10.  Thus,  for  R_  * 10 

d (3/2  u2)  _ top  u2 


which  leads  to 


^11  ^ ^ e 1 -'2 ^ (t-f) 

L 


where  Ujj  is  the  velocity  In  the  last  stage  of  decay  and  where  u is  the  transition 
velocity  determined  from 


uy  ./ 


17,  Tennekes,  11.,  and  Lumley,  J.  L.  (197  2)  A First  Course  In  Turbulence  (The 
MIT  Press,  Cambridge,  MA), 


The  time  t for  the  transition  to  take  place  can  be  determined  from  Eq.  (53) 
where  is  replaced  by  u and  determined  from  Eq.  (46) 


t ■ 3 L <4-  r~)  . 

a 


(58) 


W*  can  now  estimate  K,  At,,  from 


M 


„ t 

K°  AtM  * L / Uj  (t)  dt  + L • } ujj  <t)  dt 


(00) 


Inserting  for  Uj  and  the  expressions  in  EqB.  (54)  and  (57)  and  performing  the 
integrations  we  obtain 


K°  AtM  - 3 L*in 


✓ 3L  + u4t\  L- 

\ ~~5L — ) + T 


(61) 


where  a ■ ( 10y/3L^>.  The  second  term  on  the  right  turns  out,  under  the  substitu- 
tions below,  to  be  relatively  small. 


Inserting  typical  values  for  the  stratosphere  of  v ■ 1.  84  X It)"4  m2/Bec, 


ting 

S ■ NO^)"1'2,  N-  0.825  S*1,  L « 100  m,  and  R{  ■ 0.  15  in  Eqs.  (46),  (58),  (50), 
and  (81)  we  arrive  at 


K°  A‘m“  3,04X  l°5  m2 


(62 


using 


K°-Lut  (83 

we  find  » 308  m2/sec,  hence 

AtM  «•  1.  28  X 103  sec  (64 

||l 

which  agrees  with  Atf  * 1500  seconds.  Note,  in  passing,  that  here  far  sur- 
paaeee  the  limit  of  Eq.  (41), 

Of  course,  P+  is  related  to  Atj^  and  At^  as  we  have  seen;  however,  there  is  a 
subtle  point  to  bo  cleared  up  regarding  the  definition  of  P*.  In  the  calculations 

It  should  bu  mentioned  that,  had  it  turned  out  that  At*,  > > Atf,  the  assumption  that 
we  ignore  transport  between  the  (potentially)  turbulent  layers  could  have  been 
violated,  Once  rendered  turbulent,  they  could  have  churned  for  a very  long  time 
and  eventually  overlapped. 


associated  with  Eq,  (18)  (Ref,  I)  we  were  working  with  vertical  profiles  of  unstable 
layers  (potentially  turbulent).  In  contrast,  the  P*,  which  occurs  In  the  calculations 
by  Lilly  used  In  Eq.  (1),  refers  to  the  probability  of  encountering  turbulence  in  the 
horizontal  direction.  We,  therefore,  will  put  this  distinction  between  vertical  and 
horizontal  determination  of  P*  In  evidence  by  using  Py*  and  respectively,  The 

two  cases  of  K are  thus  rewritten  as 

c 


K. 


At.. 

p # K — -M. 


(35) 


where  the  constant  has  been  set  equal  to  1 In  Eq.  (36),  and 


K 


Py*  L2 

2 A tj. 


(2) 


We  must  re -examine  our  results  In  the  l ight  of  this  new  distinction, 

First  let  us  re-examine  the  difference  between  Atf  and  At^.  Atf  is  the  duration 
of  "potential"  rather  than  actual  turbulence  beeauae  It  represents  the  time  between 
the  event  of  the  reduction  of  R^  to  values  below  0,  25,  and  the  event  where  returns 
to  values  above  0,  25  due  to  the  effect  of  turbulent  onset  followed  by  layer  expansion 
In  other  words,  At^  is  "growth  to  turbulence  time.  " A source  of  error  in  the  model 
Is  that  sometimes  R^  dips  below  0,  25  and  then  goes  back  to  a value  above  0.  25  with- 
out an  intervening  billow  event  (see  figures  in  Browning).9  This  la  caused  by  varia- 
tions in  mean  flow,  If  this  were  taken  Into  account  it  would  show  that  a certain 
fraction  of  potentially  turbulent  layers  do  not  become  turbulent  and  hence  our  esti- 
mate of  Ke  from  Eq,  (2)  is  an  upper  limit.  One  gets  the  Impression  from  Browning's 
data  that  this  artifact  Is  not  negligible)  however,  when  one  is  more  concerned  about 
orders  of  magnitude  rather  than  factors  of  2 it  appears  perfectly  safe  to  ignore  this 
artifact.  In  contrast  AtM  Is  the  "duration"  of  the  turbulence  as  we  have  seen.  As 
mentioned  one  of  the  uncertainties  In  the  use  of  Eq,  (35)  is  the  value  to  use  for  I<L 
since  at  present  one  never  knows  A tM  or  k”  of  any  particular  In  situ  measurement. 
The  measured  value  of  may  refer  to  a much  decayed  value  for  example,  There 
are  other  uncertainties  which  will  be  discussed  elsewhere,  regarding  the  assump- 
tions In  going  from  raw  data  to  KL  (Dewan).  10  The  assumption  of  an  inertial  range  , 
for  example,  Is  not  valid,  but  the  relation  between  K.  and  t seems,  nevertheless, 
to  have  some  validity  (Panofsky  et  all. 


18.  Dewan,  E.  M,  (1079b)  Anomalous  -5/3  "Turbulence"  Spectra  iri  Stratified  FlutdB 

and  Vertical  Transport.  In  Stratosphere  1 (unpublished,  submitted  to  the 
■J,  Atm,  Sol, ). 

19,  Panofsky,  H.  A, , and  Heck,  W.  ( 1975)  Stratospheric  mixing  estimates  from  heat 

flux  measurements,  The  Natural  Stratosphere  of  1974.  CIAP  MONOGRAPH  1. 
DOT -TST -75-51  pp  8-STO  lo"8-02, — 


24 


As 


Next  we  examine  the  relations  between 
we  have  seen 


and  P,,  ;AtM>  Atf,  and  At 


V 


Pv*.  (At^/iAtjjp) 


(85) 


where,  now.  we  make  the  distinction  that  At^p  is  the  time  at  a given  altitude,  be- 
tween occurrence  of  "potentially  turbulent"  events  by  introducing  the  second  sub- 
script.  P.  In  the  case  of  P^  we  would  have  (assuming  stationarity) 


At 


'ST 


M 


'bT 


(66) 


where  AtfaT  is  the  time  between  actual  turbulent  events.  A derivation  of  Eqs.  (66) 
and  (66)  would  involve  P*  as  probabilities  for  which  At^  is  regarded  as  times 
between  "hits"  andAtf  andAtM  regarded  as  times  between  "trials"  ("frames"). 
These  times  are  the  inverse  of  probability  "frequencies,  " 

Next  we  make  the  assumption  (in  the  context  of  previously  mentioned  cautionary 
remarks) 


AtbP  “ MbT 


(67) 


and  we  turn  to  the  question  of  the  significance  of  setting  AtM  * Atf  in  the  context  of 
the  new  distinction  between  and  Py*'.  From  Eq.  (66)  and  Eq.  (66)  together 
withEq,  (67),  it  is  nowobvlous  that  this  equality  between  At^  and  Atf  is  the  same 
thing  as 

* >!» 

Py  “ PH  . (68) 

Thus  we  are  lad  to  the  question  of  whether  this  is  in  accord  with  available  experi- 
mental data,  In  Ref,  1,  P.,  -4.79  X 10"Z,  In  contrast  the  P,,*  valuea  listed  by 
- 2 v - 2 1-1 

Lilly  ranged  from  2 X 10  over  water  to  6.  2 X 10  over  high  mountains.  The 
equality  Eq.  (88)  is  therefore  not  in  conflict  with  observation)  and,  from  this  fact 
we  have  obtained  some  confirmation  that  the  theoretical  agreement  between  the 
values  of  AtM  and  Atf  (shown  above)  is  indeed  valid, 

2.6  A Lilt  of  Experimental  and  Theoretical  Queitiom  of  Intereat 
Ratted  by  the  Ke  Eitimation  Problem 

This  report  has  raised  a number  of  questions  and  it  would  be  of  convenience  if 
the  most  important  of  these  were  listed  together  in  one  place. 
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1.  AtbT  andAtbp  are  assumed  equal,  but  how  accurate  is  this  equality?  The 
first  is  the  time  between  actual  turbulent  events  at  a given  altitude  and  the  second 
is  the  time  between  events  where  < 0.  25.  It  would  seem  that  the  most  direct 
approach  to  this  would  involve  simultaneous  radar  and  radiosonde  observations  in 
the  stratosphere, 

2.  What  is  a good  value  for  Atf?  This  was  obtained  "quasi  experimentally" 
in  Ref.  1,  but  more  observations  of  the  type  in  No.  1 above  could  decrease  the  un- 
certainty of  its  value,  Atf  is  perhaps  the  most  uncertain  parameter  in  the  vertical 
stack  model  for  Kfl. 

3.  The  value  of  Py*  <and  the  associated  P (L))  is  available  for  only  one  geo- 
graphic locality  at  present.  One  needs  a more  diverse  set  of  measurements  so  that 
one  could  estimate  the  global  value  for  this  quantity.  M 

4.  A better  experimental  estimate  for  At^  of  / dt  is  needed.  In  situ 
balloon  measurements  located  within  turbulent  layer?  would  be  helpful  In  this 
regard. 

5.  Since  measurements  of  Py  ' , and  P^'  and  P (L)  and  so  on,  to  date  have 
been  made  only  under  ideal  weather  conditions,  it  would  be  useful  to  know  how  the 
estimates  are  altered  by  other  weather  conditions  such  as  storms, 

6.  Some  measurements  of  turbulence  may  have  actually  been  measurements  of 
gravity  waves.  Stewart  (1989)  has  suggested  a way  to  distinguish  between  the  two 
phenomena  and  this  should  be  investigated.  Gravity  waves  cause  no  mixing,  thus 

if  they  are  mistaken  for  turbulence,  an  overestimate  for  K would  result,  (of. 

20  0 
Dewan,  in  preparation). 

7.  All  estimates  of  K due  to  local  turbulence  to  date  have  assumed  that 

turbulent  layers  form  at  random  altitudes  and  random  times,  This  assumption  is 

so  crucial  that,  If  It  were  significantly  violated,  our  model  for  K0  would  have  to  be 

modified.*  Crane21  has  given  some  radar  data  that  could  be  Interpreted  as  evidence 

20 

that  the  assumption  is  In  fact  violated;  and  a theory  by  Dewan  In  preparation  in- 
volving trapped  gravity  waves  would,  if  supported,  imply  that  the  assumption  in 
question  is  indeed  not  valid,  This  should  be  experimentally  investigated  (perhaps 
best  by  radar  means),  and  I consider  it  to  be  the  most  important  unanswered 
question  concerning  stratospheric  turbulenoe. 


20.  Dewan,  E,  M.  (197  9a)  Stratospheric  Spectra  Resembling  Turbulence 

(unpublished,  submitted  to  Science). 

21.  Crane,  R.  K.  (1977)  Stratospheric  Turbulence  Analysis  (Air  Force  Geophysics 

Lab. ) final  report7^PGL-TH-77-6y07,  AUO  47T4'0r 

‘ The  modification  would  be  accomplished  by  a change  in  the  definitions  A tj.  and  Kj 
as  well  as  their  values. 


26 


Finally,  there  is  the  main  question  raised  by  this  report,  namely,  is  KL  large 
or  small.  In  other  words,  which  of  the  following  relations  is  the  valid  one  for  the 


estimation  of  Ke» 
Ke  * 2itf 


Ke-KLP  (1) 

If  the  previously  mentioned  experimental  evidence  is  substantiated  by  direct  measure- 
ments, then  Eq.  (1)  would  be  ruled  out. 


4.  CONCLUSION 

The  two  methods  to  estimate  Ke  from  small  scale  stratospheric  turbulence  were 
shown  to  be  two  extreme  cases  of  the  result  of  a vertical  stack  model.  In  one  case 
depends  on  turbulent  layer  thickness  and  At^  which  is  the  time  between  profiles. 

In  the  other  oase  these  geometrical  properties  are  replaced  by  a need  to  know  K^, 
the  eddy  diffusion  within  a layer. 

It  is  not  yet  known  which  of  the  two  estimates  is  correct.  To  determine  that, 
one  must  have  further  knowledge  of  KL  and  &tM  defined  in  the  text.  These  could  be 
found  by  in  situ  measurements.  Finally,  a careful  examination  of  the  vertical  stack 
model  revealed  a number  of  unanswered  questions  which  were  listed  in  order  to 
expedite  further  experimental  stratospheric  research. 
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